AO-A043  771  HONEYWELL  INC  MINNEAPOLIS  MINN  SYSTEMS  AND  RESEARCH  —ETC  F/G  12/1 

OPTIMAL  LINEAR  CONTROL  (CHARACTERIZATION  OF  MULTI-INPUT  SYSTEMS— ETC (U) 
AUS  77  C A HARVEY»  G STEIN*  J C DOYLE  N00014— 75— C— 0144 

UNCLASSIFIED  77SRC34  0NR-CR215-238-2  NL 


1 OF  1 
AD4043  77t 

m 

- 

j 

r 

i 

[ 

L 

\ 

►- 

i 

j 

END 

DAT! 

FILMED 

9 -77 

ooc 

REPORT  ONR  CR2I5-238-2 


OfFICE  ot 


RESE^^'^ 


OPTIMAL  LINEAR  CONTROL 
(CHARACTERIZATION  OF  MULTI-INPUT  SYSTEMS] 


C.A.  Harvey 
G.  Stein 
J.C.  Doyle 

Honeywell 
SYSTEMS  & RESEARCH  CENTER 

2600  RIOGNAY  PARKWAY 
MINNEAPOLIS.  MINNESOTA  55413 

CONTRACT  NQQ014-75  C-0144 
ONR  TASK  215-238 

9 AUGUST  1977 

1 December  1975  * 31  March  1977 

Approved  for  public  release;  distribution  unlimited. 


PREPARED  FOR  THE 


l-g  OFFICE  OF  NAVAL  RESEARCH  *800  N.  QUINCY  ST.«  ARLINGTON«VA«22217 


CHANGE  OF  ADDRESS 


Organizations  receiving  reports  on  the  Initial  distribution  list  should  confirm 
correct  address.  This  11st  is  located  at  the  end  of  the  report.  Any  change  of 
address  or  distribution  should  be  conveyed  to  the  Office  of  Naval  Research,  Code 
211,  Arlington,  Virginia  22217. 


DISPOSITION 

When  this  report  is  no  longer  needed.  It  may  be  transmitted  to  other  authorized 
organizations,  or  properly  destroyed.  Do  not  return  It  to  the  originator  or 
the  monitoring  office. 


DISCLAIMER 


The  findings  In  this  report  are  not  to  be  construed  as  an  official  Department 
of  Defense  or  Military  Department  position  unless  so  designated  by  other  offic- 
ial documents. 


IINCLASSIKIKI) 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (WHEN  DATA  ENTEREOI 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3,  RECIPIENT'S  CATALOG  NUMBER 


•.-TYPEOF-RtPORI/PERlOO  COVERED 


Technical  ItS 

1 ivi  - 


OR  GRANT  NUMBER(S) 


N00014-75-C-0144  / 


PROGRAM  ELEMENTT»ROJECT.T  ASK  AREA 
a.  WORK  UNIT  NUMBERS 


_9_Au£«5tS21L 

23.  NUM9CR  OF  PAGES 


Unclassified 


15*.  DECLASSIFICATION  DOWNGRADING  SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (OF  THIS  REPORT) 


Distribution  of  this  report  is  unlimited.  Reproduction  in  wi 
permitted  for  any  purpose  of  the  U.  S.  Government. 


^p^,i:ovnd  lor  ; 
^ •tiib'it'.ou 


7.  DISTRIBUTION  STATEMENT  (OF  THE  ABSTRACT  ENTERED  IN  BLOCK  20,  IF  DIFFERENT 


IB.  SUPPLEMENTARY  NOTES 


19.  KEY  WORDS  I CONTINUE  ON  REVERSE  SIDE  IF  NECESSARY  AND  IDENTIFY  BY  BLOCK  NUMBE  H| 


('tptimal  control 
Performance  index 
Synthesis  procedure 


Kigenvalues 

KiRcnvectors 

Multi -input  systems 

Mullivariable  control  systems 


•As ,y m p to ti e properties 
Attenuation  characteristics 
Automatic  control 
Control  theory 


ABSTRACT  (CONTINUE  ON  REVERSE  SIDE  IF  NECESSAR-Y  AND  IDENTIFY  BY  BLOCK  NUMBEHI 


-suits  on  characterizations  of  optimal  multi-input  systems. 

(I)  etiuivalence  of  quadratic  state  weiphtinp  to  quadratic 
es  where  m is  the  dimension  of  the  input,  (2)  characterizatitm  of 
indexes  in  terms  of  asymptotic  eipenvalues  and  ciRF-nvectors, 
lh(-  effects  of  compensators  on  hiph  fi-equency  attenuation,  and 
ocF'dure  for  constniclinp  (juadratic  pi-rformance  indexes  basi-d 
ications. 


IINGI.A.S.SIKIKI) 

SECURITY  CLASSIFICATION  OF  THIS  PAQt  (WMC  N DAI  A t N I 1 HI  HI 


REPORT  DOCUMENTATION  PAGE 

1.  report  NUMBER 

ONlicR21 5-238-2 

2.  GOVT  ACCESSION  NUMBER 

TITLE  lANO  SUBTITLE! 

d’p'riMAL  LINEAR  CONTROL 
’(CHARACTERIZATION  OF  MU 
SYSTEMS)  . 

A 

LTI-INPUT 

//' 

7.  AUTHOR(S)  / 

j C.  A.  Alarvey 

' G.  Stein 

J.  C.  Doyle 

/‘S 

9.  PERFORMING  ORGANIZATIONS  NAME/AODRESS 

Honeywell  Inc.  , Systems  and  Research  Center  ‘ 
2600  Ridpway  Parkway 
Minneapolis,  Minnesota  55413 

11.  CONTROLLING  OFFICE  NAME/ADORESS 

Department  of  the  Navy 
Office  of  Naval  Research 
.Arlington,  Virginia  22217 

J! 

14.  MONITORING  AGENCY  NAME/ADDRESS  (IF  O 

IFFERENT  FROM  CONT.  OFF.) 

FORM 
I JAN  73 


1473  EDITION  OF  1 NOV  55ISn^Oi.ETE 


'‘e  -octicn  fTj 
- “ :>Cv-f  on  pj 

n 


dSiPi3.(;:*;/;?^a^^iii[i  co^EU 


TABLE  OF  CONTENTS 


Section 


Page 


INTRODUCTION  AND  SUMMARY 

CHARACTERIZATION  OF  OPTIMAL  MULTI- 
INPUT SYSTEMS 

The  Return  Difference  Equation 

Equivalence  of  x 'Qx  to  x 'C  'Cx  in  the  Performance 
Index 

Asymptotic  Eigenvalue/ Eigenvector  Character- 
ization 

Lemma  1 

Lemma  2 

Theorem  1 

High  Frequency  Attenuation  Characteristics 

SYNTHESIS  PROCEDURE  BASED  ON  THE 
ASYMPTOTIC  CHARACTERIZATION 

Weighting  Matrix  Construction 

Ulustrative  Example 

APPENDIX  A 

REFERENCES 


rSECKDIi'D  PA'iE  tLANK-HCT 


LIST  OF  ILLUSTRATIONS 


Figure 


Page 


Closed -Loop  System  with  Compensator 
Transient  Responses 


24 

38-39 


I 


I 

1 LIST  OF  TABLES 


Table 

Page 

1 

Performance  Index  Specified  in  Terms  of 
Asymptotic  Properties 

21 

2 

F and  G Matrices 

32 

3 

C Matrix 
o 

35 

4 

Asymptotic  Eigenvalue/Vector  Characteristics 

37 

1 

V 

I 

^ L.  . 


SECTION  I 


INTRODUCTION  AND  SUMMARY 


This  research  program's  goal  is  the  advancement  of  optimal  linear 
control  design  technology  to  the  point  where  it  yields  operational  control 
laws  that  meet  conventional  design  specifications.  The  results  for  single 
input  systems  are  reported  in  [1],  This  report  presents  results  on 
the  characterization  of  optimal  multi -input  systems. 

The  major  results  established  during  this  phase  of  the  program  were: 

• Equivalence  of  quadratic  state  weighting  to  quadratic  weighting 
of  m responses  where  m is  the  dimension  of  the  input, 

• Characterization  of  quadratic  performance  indices  in  terms 
of  asymptotic  eigenvalues  and  eigenvectors, 

• Characterization  of  the  effect  of  compensators  on  high 
frequency  attenuation  properties,  and 

• Development  of  a procedure  for  constructing  quadratic 
performance  indices  based  on  modal  design  specifications. 

The  equivalence  result  provides  a direct  generalization  of  the  concept  of 
model  response  weighting^  for  single  input  systems  and  hence  serves  as 
an  aid  in  the  interpretation  of  the  asymptotic  properties  of  optimal  systems. 
This  result  also  provides  a general  form  for  performance  indices  which 

have  fewer  parameters  than  the  state  form  and  hence  possess  less 

. . [2J 

redundancy. 
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The  asymptotic  eigenvalue/ eigenvector  characterization  provides 
important  insight  for  the  weight  selection  process  of  linear  optimal 
design.  This  characterization  provides  a complete  description  of  the 
quadratic  performance  index  in  terms  of  nm  independent  optimization 
parameters  where  n is  the  dimension  of  the  state  vector. 

High  frequency  attenuation  properties  are  especially  important  for 
operational  control  systems.  Optimal  systems  without  compensators 
are  shown  to  have  first  order  attenuation  characteristics.  This  may  be 
insufficient  for  many  applications.  The  characterization  of  compensator 
effects  on  attenuation  provides  a guide  for  selection  of  the  proper  control 
system  structure  to  meet  design  specifications. 

The  above  results  are  derived  in  Section  II.  Section  III  presents  the 
development  of  the  procedure  for  constructing  quadratic  performance 
indices  from  design  specifications.  Application  of  this  procedure  to  an 
illustrative  example  is  also  described  in  Section  III. 

Since  transmission  zeros  of  multivariable  systems  play  a key  role  in  the 
asymptotic  characterization.  Appendix  A summarizes  the  important 
concepts  of  multivariable  system  zeros. 


SECTION  II 


CHARACTERIZATION  OF  OPTIMAL  MULTI-INPUT  SYSTEMS 


We  consider  the  linear  controllable  system 

X = Fx  f Gu.  x(0)  = x^  (1) 

where  x is  an  n-vector  and  u is  an  m -vector*  with  m generally  greater 
than  one  and  less  than  or  equal  to  n.  We  are  interested  in  basic 
characteristics  of  linear  control  laws  that  are  optimal  with  respect  to  a 
quadratic  performance  index  of  the  form 

J ■ r ” [x*  Qx  + u'  Rujdt,  Q > 0,  R > 0 (2) 

Jo  — 

with  (Q,  F)  observable.  It  is  well  known  that  the  optimal  control  is  given 

by 

u = Kx  (3) 

where  the  gain  matrix  K satisfies 

K = -R'^G'  P (4) 

with  P being  the  positive  definite  symmetric  solution  of  the  Riccati 
equation 

PF+  F'  P - PGR'^G'  P + Q = 0 (5) 


*The  components  of  u are  assumed  to  be  independent,  i.  e.  Rank  (G)  - m. 
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We  will  first  present  the  known  characterization  in  terms  of  the  return 
difference.  Then  we  present  a new  result  on  the  equivalence  of  quadratic 
performance  indices  of  the  form  (2)  to  the  form 


J =f  (r^  r " u*  Ru)dt  (6) 

J o 

where 

r = Cx  (7) 

is  an  m -vector.  Optimal  systems  are  then  characterized  in  terms  of 
asymptotic  eigenvalues  and  eigenvectors.  A unique  relationship  between 
these  asymptotic  properties  and  the  weighting  matrices  in  the  performance 
index  is  then  established.  Finally  optimal  systems  are  characterized  in 
terms  of  their  high-frequency  attenuation  properties. 

THE  RETURN  DIFFERENCE  EQUATION 

The  Riccati  equation  (5)  may  be  rewritten  using  (4)  as 

P(sl-F)  + (-sI-F')P  + K'RK  = Q (8) 

Multiplying  by  G'  (-sI-F')  ^ on  the  left  and  (sI-F)  on  the  right,  gives 
G'(-sI-F')’^PG  + G'P(sI-F)“^G  + G'(-sI-F')'^K'RK(sI-F)’^G 
= G'(-sI-F')'^  Q (sI-F)'^G  (9) 


J 


1 


I 


Since  PC  = -K^  R,  adding  R to  both  sides,  and  rearranging  yields* 
lI-K(-sI-F)‘^Gj'  R[I-K(sI-P)'^Gj 

= R f G'(-sI-F  )"^  Q(sI-F)'^G  (10) 

where  the  return  difference  matrix  is  the  mxm  matrix 

T(s)  = I-K(sI-F)'V.  (11) 

For  s = ju),  the  last  term  of  the  sum  on  the  right  hand  side  of  (10)  is  non- 
negative  for  all  real  Thus,  optimal  systems  may  be  characterized  by 
the  inequality 

lT(-jo))]'RT(ju))>  R for  aU  real  uj  (12) 

EQUIVALENCE  OF  x'Qx  TO  x'C'Cx  IN  THE  PERFORMANCE  INDEX 

Assume  that  we  are  given  the  system  (1)  and  performance  index  (2).  Let  us 
first  make  the  change  of  control  variables. 

V = /Ru 

where  /R  is  the  positive  definite  symmetric  square  root  of  R. 
and  (2)  may  be  rewritten  as 

A 

X = Fx  t-  Gv,  x(0)  = x^ 

J = r "ix'Qx  f v'vjdt 
o 

''  -1 

where  G = G(^R)  . The  optimal  control  is 

*I  is  used  in  equation  (10)  and  subsequent  equations  to  denote  the  identity 
matrices  of  appropriate  dimensions. 


(13) 

Then  (1) 

(14) 

(15) 


h 


V = Kx  = (/R)Kx 

The  return  difference  matrix  for  the  system  (14)  with  control  (16)  is 

T(s)  = I - K(sI-F)'^G  (17) 

and  the  return  difference  equation  is 

lT(-s)J'  T(s)  - 1 = G'(-sI-F')  ^ Q(sl-F)  ^G  (18) 

For  s = iu),  the  right  hand  side  of  (18)  is  non-negative  for  all  real  u). 
According  to  Molinari  [2],  Popov  [3J  has  shown  that  iT(-jou) J'T(ja)  - I 2;  0 
implies  the  existence  of  an  mxm  matrix  C such  that 

lT(-s)J'T(s)-I  = G'(-sI-F')’^  C'C(sI-F)G  (19) 

For  such  a matrix  C,  let 

V = (^9) 

be  the  optimal  control  for  the  system  (14)  and  the  performance  index 

T = r ” (x'C  tx  + v'v)dt  (31) 

o 

The  return  difference  equation  for  the  system  (14)  with  the  control  given 
by  (20)  is 

[T(-s)J'f(s)-I  = G'(-sI-F')'^  C'C(sI-F)'^G  (22) 

where 

T(s)  = I-?^(sI-F)‘^G.  (33) 
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Thus,  from  (17),  (19),  (22)  and  (23),  we  have 

[T(-s)J'T(s)  = II  - k(-sI-F)'^Gj 'II  - k(sI-F)'^Gj 

= II  - T^(-sI-F)'^Gj'II  - R(sI-F)'^GJ  (24) 

A A 

The  factorization  of  lT(-s)J  T(s)  was  shown  by  Youla  1 4,  Theorem  2j  to  be 
unique  to  within  a constant  unitary  matrix,  which  implies 

I - k(sI-F)'  G = MU  - K(sl-P)  G]  (25) 

where  M is  constant  and  M'M  = I.  Equation  (25)  may  be  rewritten  as 

I-M  = (k-MK)  (sI-F)'^G  (26) 

The  left  hand  side  of  (26)  is  independent  of  s,  and  the  limit  of  the  right 
hand  side  of  (26)  as  s tends  to  infinity  is 

lim  (k-MK)  6/s  = 0 (27) 

S-^co 

which  implies  M = I,  Thus, 

0 = (k  - k) (sI-F)'^G  (28) 

which  is  an  identity  in  s.  For  |s|  large,  (sI-F)  ^ may  be  written  as 

03 

(sI-F)'^  = s'^  E (Fs'^^  (29) 

j = o 

Therefore, 

(k  - k)F^G  = 0,  j = 0.  1,  2.  ...  (30) 


Equation  (30)  together  with  the  controllability  of  (F,6)  implies  tliat  K = K. 
This  establishes  that  the  performance  indices  (15)  and  (21)  are  equivalent 
in  the  sense  that  they  yield  the  same  optimal  controller.  Writing  these 


■j 
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performance  indices  in  terms  of  the  original  control  vector,  u,  gives  the 
equivalence  of  the  performance  index  (2)  to  (6)  for  the  system  (1)  with  r 
defined  by  (7)  for  some  mxn  matrix  C. 


ASYMPTOTIC  EIGENVALUE/ EIGENVECTOR  CHARACTERIZATION 

For  this  discussion  we  consider  the  system  (1)  and  the  quadratic  perfor- 
mance index  in  the  form 

J = r”  (r'r+  pu'Ru)dt  (31) 

where  r is  given  by  (7)  and  p is  a positive  real  scalar  parameter.  We  are 
interested  in  the  asymptotic  properties  of  the  optimal  closed -loop  system 
as  the  parameter  p tends  to  zero.  A complete  characterization  will  be 
given  for  the  case  when  the  pair  (C,  F)  is  observable  and  the  matrix  CG 
has  rank  m. 

The  optimal  control  is  u = Kx  where  the  feedback  gain  matrix 
K = -p'^  r'^G'P 

and  P = P^  > 0 is  the  solution  of  the  Riccati  equation 
PF  + F'P  - p'^PGr'^G'P  + C'C  = 0 

The  return  difference  equation  for  the  optimal  system  is 

lT(-s)J'pR  T(s)  = pR  f G'(-sI-P')’^  C'C(sI-F)“^G  (34) 

where  T(s)  is  the  return  difference  matrix  (11). 


is 

(32) 


(33) 


.J 
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The  determinant  of  the  return  difference  matrix  is 

Taking  the  determinant  of  (34)  and  rearranging  yields 

<t  (s)  0 (-S)  = 0 (s)  0 (-S)  I I + — r'^  H'(-s)  C'CH(s)  1 (36) 

C C O O ' p ' 

where 

0^(s)  = Isl-F|,  0^(s)  = [sI-P-GK  1 

and  H(s)  = (sI-F)  ^G.  The  polynomials,  and  0^(s),  are  the 

open -loop  and  closed -loop  characteristic  polynomials  respectively. 
Equation  (36)  may  be  analyzed  to  determine  the  asymptotic  nature  of  the 
closed -loop  eigenvalues  as  p tends  to  zero. 

The  right  hand  side  of  (36)  is  a polynomial  in  p ^ which  may  be  written  as 

4/ is)  (a)  ? a.(s)p"^  (36) 

o o . j 

J = o 

where  a = 1 and  for  j = 1,  2,  ....  m, 
o 

a^(s)  = ordered  principal  minors  of  R ^H'(-s)C'CH(s)J  (37) 

2 

We  note  that  each  is  a ratio  of  polynomials  in  s with  the  denominator 

equal  to  \I/  (-s)  \j/  (s).  Thus 
o o 

4/  (s)  4^  ("s)  a.(s)  = n.(s^)  (38) 

O ^o  3 J 

with  the  maximum  degree  of  n.  being  n-j. 
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Kwakernaak  describes  the  asymptotic  properties  of  zeros  of  general 
polynomials  of  the  form  (36)  as  p tends  to  zero  in  reference  5.  Let  us 
here  consider  the  special  case  in  which  the  degree  of  n^  is  n-m  and  the 

rank  of  CG  is  m.  In  this  case,  as  p tends  to  zero,  n-m  roots  of  (s) 

/ 2 ^ 
approach  the  roots  of  ) which  have  negative  real  parts.  These  are 

the  zeros  or  the  left  half  plane  mirror  images  of  the  zeros  of  the 

determinant: 

lC(sI-F)'^G|  (39) 

The  zeros  of  (39)  are  the  transmission  zeros  associated  with  the  input- 
output  system  described  by  (1)  and  (7).  Appendix  A summarizes  recent 
research  on  system  zeros.  The  remaining  m closed -loop  eigenvalues 
approach  infinity  in  multiple  Butterworth  patterns.  To  examine  the 
nature  of  these  eigenvalues,  we  consider  0^(-s)  0^(s)  = 0 for  |s  | 
sufficiently  large  so  that  0^(-s)  0^(-s)  ^ 0.  Then  from  (36)  we  are 
interested  in  the  equation 

0 = 1 1 + p"^r'^H'(-s)  C tH(s)| 

= I I + p'^R'^hVs)  C*CH(s)r’^  I (40) 

For  |s  I large  (40)  may  be  written  as: 

ll  - -s'^  (D  +0(|sf^))l  = 0 (41) 

P 

where  D = r'^G'C'CGR*^ 

-2  I I -2 

and  0(|sl  ) indicates  terms  which  tend  to  zero  as  |s|  -♦  ®. 
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Thus 

= - X(D)+  0(|s  r^)  (42) 

P 

where  x(Matrix)  denotes  the  eigenvalues  of  the  matrix,  and  in  our  case 

_2 

D > 0 so  that  X (D)  are  real  and  positive.  Hence  x[D  + 0(|s|  )]  = 

_2 

X (D)[l  + 0(|s|  )]  so  that  as  p - 0,  m poles  ap^jroach  infinity  along  the 

' ' -k 

real  axis  at  the  rate  of  p ® . 

Kwakernaak  [5j  points  out  that  if  the  determinant  (39)  has  less  than  n-m 
zeros  or  if  the  rank  of  CG  is  less  than  m then  the  eigenvalues  which  go  to 
infinity  may  group  into  several  Butterworth  patterns  with  different  growth 
rates. 

It  is  somewhat  disconcerting  that  when  C is  of  rank  m at  most  n-m 
asymptotic  eigenvalue  locations  may  be  specified  in  the  finite  plane.  It 
is  possible  to  have  more  than  n-m  finite  asymptotic  eigenvalue  locations 
by  choosing  C to  have  rank  less  than  m.  Such  an  input-output  system 
is  called  degenerate  (see  Appendix  A).  If  C is  chosen  to  be  a qxn  matrix 
where  q < m,  then  equation  (36)  gives 

0 (s)0  (-S)  = 0 (s)0  ( -s)  II  + — H'(-s)C'CH(s)  1 
c ^c  o o ' m p ' 

= 0^(s)0^(-s)  1 Iq  + ^CH(s)  r'^  H'(-s)C'1 

= 0 (s)0  (-s)  f I + "*  + p"'’  CH(s)r’^  n'(-s)C'|  -i 
° ° (43) 
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The  elements  of  CH(s)R  ^ are  ratios  of  polynomials,  the  denominators 

being  of  degree  n and  the  numerators  of  degree  less  than  or  equal  to  n-1. 

The  determinant,  | CH(s)R  ^H'(-s)C  | is  a ratio  of  polynomials,  the 

denominator  being  of  degree  2qn  and  the  numerator  being  of  degree  at 

most  2q(n-l).  The  degree  of  t//^(s)i//^(-s)  is  2n.  Cancellation  of  the 

denominator  with  \1/  (s)i//  (-s)  and  common  factors  in  the  numerator  takes 
o o 

place  in  (43).  Thus  the  coefficient  of  p in  (43)  is  a polynomial  of 
degree  at  most: 

2q(n-l)  + 2n  - 2qn  = 2(n-q)  (44) 


The  matrix  C may  be  chosen  to  achieve  the  maximal  degree,  2(n-q).  In 
this  case  there  are  n-q  finite  asymptotic  pole  locations.  It  appears  that 
this  may  be  of  some  advantage  as  it  allows  more  poles  to  be  specified. 

If  a vector  C were  used,  then  q=l  and  it  is  possible  to  have  n-1  finite 
asymptotic  poles.  Unfortunately,  ingeneral  whenq<m,  all  n-q  poles 
may  not  be  specified  arbitrarily.  Consider  the  following  example 
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0 

1 

0 

0 

0 

0 

0 

0 

0 

1 

- 

. 

_ 

R 


C 


tc, 


^2  ^3  ^^4^ 


12 


r 


Then 


8.  -1  8_ 
s +p  s C 


-2—0 
s s 


i-2  ^ 

s s 


8 -1  8_ 
= s + p s C 


1 . -1 


3 0 0 


-3  -2  0 0 


^^4  -^3 

s s 


4,  4 . -1  2 


-3  -2 

s s 


{s’+  p " [(C^'  + Cg")  -(c^  )s^l 


Thus  while  three  poles  will  be  finite  as  p -*  0,  two  will  be  at  zero 
independent  of  the  choice  of  C.  In  this  example  the  pair  (C,  F)  is  not 
observable,  but  similar  examples  can  be  constructed  with  an  observable 


No  general  characterization  of  the  nature  of  the  pole  restriction  has  been 
found.  Apparently  the  only  general  statement  we  can  make  is  that  for  a 
qxn  matrix  C,  where  q<m,  the  n-q  finite  poles  cannot  be  arbitrarily 
specified. 
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For  an  mxn  matrix  C it  is  well  known  that  the  zeros  of  (39)  may  be  selected 
arbitrarily  by  appropriately  choosing  C.  There  are  a multitude  of  methods 
for  selecting  C to  achieve  desired  zeros.  The  difficulty  is  that  the  choice  of 
C is  not  unique  since  only  n-m  parameters  are  needed  to  specify  the  n-m 
zeros.  This  difficulty  with  nonuniqueness  will  be  removed  by  considering 
the  freedom  to  specify  eigenvectors  corresponding  to  these  asymptotic 
eigenvalues.  The  characterization  will  be  completed  with  consideration  of 
the  m eigenvalues  that  approach  infinity  and  their  associated  eigenvectors. 

We  begin  the  discussion  of  eigenvectors  with  the  following  result  established 
by  Moore  [6]: 

Lemma  1 


The  vectors  x.  and  associated  distinct  complex  numbers  s.  are 
1 1 

closed-loop  eigenvectors  and  eigenvalues,  respectively,  of  the 
system  (1)  with  feedback  (3)  if  and  only  if  there  exist  m -vectors 
such  that 

(s.I-F)'V.  M-  = X.  (45) 

1 1 

= Kx^  i = 1,  2,  ...,  n (4(i) 

This  lemma  shows  that  the  eigenvectors  achievable  with  linear 
state  variable  feedback  are  confined  to  m -dimensional  subspaces 
determined  by  their  associated  eigenvalues.  This  eigenvector 
adjustment  plus  arbitrary  pole  placement  is  precisely  the  freedom 
offered  by  linear  state  variable  feedback. 

The  m-vectors,  up  rnay  be  characterized  with  respect  to  the  return 
difference  as  follows: 


14 


Lemma  2 


V 

i 

I 

I 


) 

I 


I 


The  m-vectors  are  identical  (except  possibly  for  magnitude) 
to  the  nonzero  vectors  determined  by 

T(s^)  = 0,  i = 1,  2,  . . . , n (47) 

To  show  this,  note  that  the  ^i,^'s  from  (45)  and  (46)  satisfy  (47),  i.  e. 

T(s.)  p,.  = ti,.  - K(s.l-F)’^G  u. 

1 1 1 

= u-  * Kx- 
1 

= 0 (48) 


Moreover,  solutions  of  (47)  are  unique  except  for  magnitude  as 
long  as  the  s^'s  are  distinct.  For  this  case,  the  s^'s,  interpreted 
as  transmission  zeros  of  T(s),  have  unit  geometric  and  algebraic 
multiplicity  and,  hence,  the  rank  deficiency  of  T(s)  at  s^  is  one. 

It  then  follows  that  with  appropriate  normalization 

i = 1,  2,  . . . . n (49) 


For  the  case  in  which  the  rank  of  CG  is  m.  we  may  sumniai’izc  the 
asymptotic  eigenvalue/eigenvector  properties  in  the  form  of  a theorem. 
With  the  rank  of  CG  being  full  and  hence  the  rank  of  the  mxn  matrix 
C being  full  there  is  no  loss  in  generality  to  normalize  C to  be  of  the 
form 


c = w c 

o o 

where  is  a nonsingular  mxm  matrix  and 


C 

o 


I J 


m 


(50) 


(51) 
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Such  a normalization  maybe  achieved  by  a state  transformation  which  reorders 
the  state  variables.  With  this  notation  we  may  state  the  following  theorem. 

Theorem  1 

Let  r = Cx  with  C satisfying  the  conditions 

i)  the  rank  of  CG  is  m 

ii)  C is  in  the  form  described  by  equations  ) and  (51) 

iii)  the  zeros  of  |C(sI-F)  ^G|  are  distinct  and  have  negative 
real  parts. 


Then,  the  optimally  controlled  system  (1)  with  respect  to  the  performance 
index  (31)  has: 

n-m  asymptotically  finite  eigenvalues 
s^°,  i = 1.  2,  n-m  (52) 

and  m asymptotically  infinite  eigenvalues 


s.  p i = 1,  2,  ....  m (53) 

where  the  s^°  have  negative  real  parts  and  if  for  some  i,  is  not  real 
there  is  a s.°  which  is  the  complex  conjugate  of  s°  and  the  are  real  and 
negative.  The  s^°'s  are  distinct  and  have  associated  eigenvectors 

x.°  = (s.°I-F)‘^G  v.°  (54) 

where  the  s.°  and  v.°  are  defined  by 

C(s.°I-F)‘^G  v.°  =0.  y 0 (55) 


16 


The  eigenvectors  associated  with  the  asymptotically  infinite  eigenvalues 


X.  = G V. 

1 1 

CD  CO 

where  the  s^  and  are  defined  by 

R = (N')'^  S'^n"^ 

W'W  = [(C  GN)'J'Mc  GN)'^ 
o o o o 


A r CO  00  00  - 

N ^ [v.  v„  ...  V ] 

1 J m 

S - diag  (s,  s„  ...  s ) 
I ^ m 


This  theorem  combines  known  results  from  the  literature 


[5.7] 


concerning 


the  asymptotic  behavior  of  eigenvalues  with  new  results  for  the  behavior 
of  the  corresponding  eigenvectors.  The  following  proof  emphasizes  the 
new  results. 

Proof  of  Theorem: 

From  Lemma  2,  the  optimal  controller's  eigenvalues  and  eigenvectors 
are  characterized  by  (47)  which  may  be  premultiplied  by  pT  (-s^)R  yield- 
ing 

pT'(-s.)R  T(s.)  V.  = 0 (61) 

^ 1 11 

Then  from  the  return  difference  equation  (34),  we  have 

[pR  + H^(-s.)  C'CH(s.)Jv.  = 0 (6  2) 
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We  are  interested  in  the  asymptotic  behavior  of  the  pairs  (s^,  v.)  as  p 
tends  to  zero.  In  [7J  it  is  shown  that  if  the  rank  of  CG  is  m,  then  there 
will  be  2(n-m)  such  pairs  for  which  the  s^  values  remain  finite.  The  s^ 
values  will  be  symmetric  about  the  imaginary  axis.  Those  in  the  left 
half -plane  will  be  the  eigenvalues  of  the  stable  closed -loop  system. 

Taking  the  limit  of  ((>2)  as  p approaches  zero  yields 

H'(-s.)  C'CH(s.)v.  = U (fi3) 

1 11 

Condition  iii)  and  the  symmetry  condition  on  the  s^  values  imply  that 

[II^(-s.)C'|  A 0 for  s.  with  negative  real  parts.  Hence  from  ((i3)  we  have 

for  the  stable  s., 

1 

C H(sJv.  = 0,  V.  =/  0,  i = 1,  2,  . . . , n-m  (ti4) 

which  establishes  equation  (55),  and  equation  (54)  follows  from  Lemma  1. 

I 5 J 

It  is  known  that  the  remaining  m closed -loop  eigenvalues  approach 
infinity  in  m first-order  Hutterworlh  patterns.  Their  corresponding 
eigenvectors  can  be  determined  by  using  a power  series  expansion  in  s. 
of  equation  (62),  i.  e. 

0 =(pR  - (s.)'^[G'(I-F's.'^  f ...)C'JC(H-F  s.'^  t . . . )G]  v. 

= p{R  - (p^  s.)'^[(CG)tG  + 0(s.'^)J]  V. 

= p[R  - (a.)""^  I(CG)'CG  f 0(pa.'^J}  V.  (65) 

where  a.  = s.  p‘  . As  p tends  to  zero,  a.  approaches  s.  and  the 
asymptotically  infinite  pairs  satisfy 

0 = [R  - (s  ")*^  (CG)'CGjv.",  V."  y 0 (66) 

O 11 
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pmtm 


5 


From  this  equation  and  condition  i)  we  can  derive: 


0 = {[(CG)'J'^  R(CG)'"  - (sp'^I}CG  V.”  v."  A 0.  i = 1.  2, 

(67) 


,-l 


m 


This  equation  implies  that  the  symmetric  matrix  [(CG)^J  R(CG)  has 


uo  - 2 

eigenvalues  (s^”)  and  orthogonal  eigenvectors  CG  If  we  normalize 


the  by  the  condition 


jCG 


= 1,  i = 1,  2,  . . . , m 


we  have 


(CGN)'CGN  = I 


(68) 


(69) 


and 


-2 


R = (CG)'CGN  S " (CGN)'CG  (70) 

where  the  matrices  N and  S are  defined  in  (59)  and  (60).  From  condition 
ii),  equations  (69)  amd  (70)  may  be  rewritten  as 

W'  W = {(C  GN)'J'^  (C  GN)'^ 
o o o o 


and 


R = (C  G)'  W 'W  C GNS'^(C  GN)'W  'W  C G 
o ooo  o ooo 


= (C  G)'[(C  GN)'j'^s'^(C  GN)  G 
o o o o 

= (N')’^  S‘^n“^ 

establishing  equations  (58)  and  (57).  Equation  (56)  which  characterizes 
the  corresponding  eigenvectors  follows  directly  from  Lemma  1 for  s^  = 

p ^ tending  to  infinity  in  magnitude  as  p tends  to  zero.  This  completes 
the  proof  of  Theorem  1. 


,-L 
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This  theorein  providos  important  now  insights  for  the  weight  selection 
process  in  linear  optimal  design. 


We  note  first  that  the  asymptotic  properties  provide  a complete  character- 
ization of  the  performance  index  (31)  in  terms  of  nm  free  parameters.  The 
finite  modes  are  characterized  by  (n-m)  eigenvalues  (s^°)  plus  (m-1)  (n-m) 

parameters  for  their  associated  normalized  directions  (v  °).  We  will 

1 

show  in  Section  III  that  these  provide  a unique  definition  of  via  equation 
(55). 

Similarly,  the  asymptotically  infinite  modes  are  characterized  by  (m-1) 
parameters  for  the  normalized  eigenvalues  (s.“)  and  (m-1)  (m)  for  their 
associated  normalized  directions  (v^  ).  These  provide  unique  definitions 
of  and  W^  via  equations  (57)  and  (58).  Together  with  the  scalar  p, 
these  parameters  specify  the  criterion  completely.  These  observations 
are  summarized  in  Table  1.  Note  that  the  scalar  p is  the  only  remaining 
"trade-off"  parameter  in  the  specification.  It  indexes  the  optimal  eigen- 
value and  eigenvector  loci  implicit  in  C , R and  W , thus  trading  off 
control  energy  against  the  degree  to  which  the  asymptotic  properties  are 
achieved. 

We  note  next  that  the  asymptotically  finite  eigenvalues  (s.°)  correspond  to 

1 8 J ^ 

transmission  zeros  of  the  square  response  transfer  matrix  CH(s)  and 

that  the  vectors  (x^°,  v^°)  are  "state-  and  control-zero  directions  [9J" 

associated  with  these  transmission  zeros.  Moreover,  for  each  fixed 

finite  s^°,  the  eigenvectors  achievable  with  linear  optimal  design  correspond 

precisely  to  those  achievable  with  arbitrary  linear  state  feedback  (Lemma 

1).  The  existence  of  zeros  and  zero  directions  in  CH(s)  are  algebraic 


I 


Table  1.  Performance  Index  Specified  in  Terms  of 
Asymptotic  Properties 


structural  properties  of  the  performance  index  which  must  be  "built  into" 
the  criterion  if  we  wish  to  achieve  desirable  dominant  modes  and  mode 
distributions  in  the  final  closed -loop  system. 

As  a final  observation,  we  note  that  Theorem  1 provides  for  the  first  time 
an  explicit  way  to  adjust  the  behavior  of  asymptotically  inifinite  modes. 

'fhese  modes  have  particular  significance  to  designers  because  ttiey 
govern  the  characteristics  of  control  actuation  devices.  Note  that  the 
eigenvectors  of  asymptotically  infinite  modes  are  in  the  range  spact?  of  G. 
Hence,  if  the  state  model  includes  actuator  dynamics  with  G nonzero  only 
for  actuator  rows,  the  asymptotically  infinite  modes  are  actuator  modes. 
According  to  equations  (53)  and  (56)  - (60),  these  modes  have  bandwidth 
ratios  and  cross  couplings  specified  entirely  by  computable  matrices 
and  W^.  In  the  past,  designers  were  forced  to  struggle  iteratively  to 
achieve  bandwidth  ratios  and  couplings  consistent  with  hardware  constraints. 

HIGH-FREQUENCY  ATTENUATION  CHARACTERISTICS 

We  will  first  show  that  optimal  multi-input  systems  have  the  property  that 
the  return  ratio  matrix  at  the  input  possesses  a first  order  attenuation 
characteristic  for  large  frequencies.  Then  we  show  that  including  a compen- 
sator increases  the  order  of  attenuation  at  the  actuator. 

For  the  system  (1)  with  control  being  optimal  with  respect  to  (2),  i.e. 
the  control  u is  determined  by  (3),  (4),  and  (5),  the  return  ratio  matrix  at 
the  input  a is: 

A^(s)  = K(sI-F)'^G 


(71) 


Expanding  the  right  hand  side  of  (71)  in  a power  series  in  s ^ for  |s| 
large  yields 

00 

A (s)  = E s’  F-’g  (72) 

“ j=o 

The  matrix  KG  = -R  ^ G'PG  is  nonsingular  since  R and  G'PG  are  positive 
definite.  Thus,  we  may  write  (72)  as 

A^(s)  = s'^  KG  [I  + 0(s'*)J  (73) 

for  js|  large.  This  demonstrates  that  A^(s)  possesses  the  first  order 
attenuation  characteristic. 


Now  let  us  consider  the  case  in  which  a compensator  is  included.  For  this 
purpose  let  us  suppose  the  vector  u consists  of  two  subvectors  u^  and  u^ 
and  that  the  compensator  is  inserted  in  the  U2  path.  The  block  diagram 
of  such  a system  is  shown  in  Figure  1,  The  differential  equation 
description  for  the  system  is 

i = Fx  + GjU^  + G^u^ 


X = F X + G V 
c c c c 

"2  ’ '“‘o 


<74) 


The  performance  index  for  this  system  is 


■^11 

( 

0*^ 

X 

+ lu^.v'J 

■ «11 

^2' 

.<^'12 

<^22j 

X 

c 

«22. 

V 

}dt 
(75) 


I 


23 


Figure  1.  Closed-Loop  System  with  Compensator 


1 


The  remaining  terms  in  may  be  written  for  large  (s|  as 

<V2i  = iiW-i.y-‘c^U-‘K^,Gj 


Equations  (85)  and  (88)  show  that  the  compensator's  attenuation 
characteristic  is  " added"  to  the  first  order  (at  least)  characteristics  of 


the  u^  subvector. 


SECTION  III 


SYNTHESIS  PROCEDURE  BASED  ON 
THE  ASYMPTOTIC  CHARACTERIZA  ITON 


In  ttiis  section  we  will  describe  a procedure  for  constructing  the  quadratic 
weighting  matrices  of  the  performance  index  from  desired  sp(,*cifications 
of  closed-loop  eigenvalues  and  eigenvectors.  As  an  illustrative  example 
the  results  obtained  with  this  procedure  for  the  lateral  axis  controller 
design  for  the  F-4  will  be  described. 

WEIGHTING  MATRIX  CONSTRUCTION 

A key  feature  of  the  asymptotic  interpretations  above  is  the  structural 
requirement  that  C(sl-F)  ^G  must  possess  transmission  zeros  and  zero 
directions  corresponding  to  the  pairs.  Tliis  section  provides  a 

construction  procedure  to  "build"  such  response  transfer  matrices  based 
on  desired  closed-loop  specifications. 

We  begin  with  the  assumption  that  (without  loss  of  generalil>  ) llu'  matrix 
G has  been  transformed  to  the  form 


with  G.^  an  mxm  nonsingular  matrix. 


Further,  we  assume  that  for  desired  closed-loop  eigenvalues,  s^  , the 
corresponding  vectors  x.  which  describe  the  desired  mode  distributions 
in  state  space^  can  be  specified  by  the  designer,  i.  e.,  the  desired  modes 
should  look  like 

x^(t)  = e i ; i = 1,  2,  ...,  n-m  (88) 

* 

It  will  generally  be  true  that  only  a few  of  the  components  in  x^  are 

actually  specified.  The  rest  can  be  arbitrary.  To  account  for  this,  we 

♦ 

reorder  and  partition  x^  as  follows: 


♦ «i  IV 


where  v is  the  specified  subvector,  v's  are  the  unspecified  components, 
and  {.. .}  denotes  the  reordering  operation. 


In  accordance  with  Lemma  1 and  Theorem  1,  the  desired  mode  (89)  may 
not  belong  to  the  set  of  eigenvectors  achievable  by  optimal  linear  design; 


o , R.  , 

V;  r/  O,  , 1 O L O 

= f(s.  I-F)  G]  V.  = ---V. 

'•1  Ml 

V 


The  designer  is  free  to  select  so  as  to  best  approximate  with  . 
One  way  to  do  this  is  by  orthogonal  projections. 


v.°=  (r/u'^r/y.* 


A known  mode  distribution  in  some  output  space  y = Cx  with  C non- 
singular  would  obviously  servo  as  well. 


29 


although  other  methods  may  be  more  suitable  in  particular  applications. 


We  will  assume  here  that  ; i = 1.  2,  ....  n-m]  have  been  determined 


by  some  procedure  and  result,  via  equation  (54),  in  the  following  eigen- 
vectors 


X = lx  °x  °...  x°  J 
1 I n-m 


(92) 


where  care  was  taken  to  assure  that  no  linear  combination  of  the  x.  's  lie 


in  the  range  space  of  G.  We  then  define  the  projection  matrix 


P = I -X  (X'X)'^  X' 
n 


(93) 


and  note  that 

PG  p,  /=  0 if  p.  /=  0 


(94) 


Let  us  partition  the  matrix  P in  the  form 
P = 


p. . 

P.  „ 

11 

12 

p ' 

P .o 

12 

22 

(95) 


with  P^2  being  an  mxm  matrix.  Conditions  (87)  and  (94)  and  the  non- 


singularity of  G imply  that  P„„  is  nonsingular.  Thus,  C may  be  defined 
^ A w o 

as 


[oIp 


-1 


22 


JP  = IP 


■1, 


/ I 


22 


J 


12i  m 


(9  6) 


To  verify  that  has  the  desired  structural  property,  we  note  that 


C (s.°I-F)‘^G  v.° 

0 1 1 


[0;p^2'^JP  (s,°I-F)'V, 


1 

o 


= lo;P22  IP»1 

= 0,  i = 1,  2,  . . . , n-m 


(97) 
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and  that  for  sufficiently  large  |s|,  the  rank  of  C^(sl-F)  G is  m.  Thus,  there 

is  a rank  reduction  at  sP  which  means  sP  is  a transmission  zero,  and  hence 

1 1 

an  asymptotically  finite  eigenvalue  with  corresponding  eigenvector  x,  . 


ILLUSTRATIVE  EXAMPLE 

For  purposes  of  illustration  we  consider  the  inner  loop  lateral  axis  design 
problem  for  the  F-4  fighter  aircraftat  alow  dynamic  pressure  flight  condi- 
tion taken  from  reference  flO],  The  dynamics  are 

X = Fx  + Gu 


stability  axis  roll  rate 
stability  axis  yaw  rate 
angle  of  sideslip 
bank  angle 


rudder  deflection 


aileron  deflection 


rudder  command 


aileron  command 


Matrices  F and  G are  given  in  Table  2.  Note  that  tlicy  include  first  order 
actuator  dynamics  with  10  and  5 rad/sec  bandwidths.  These  have  been 
made  deliberately  slower  than  the  hydraulic  actuators  available  on  the  1'  4 


:n 


(approximately  20  and  10  rad/sec)  in  order  to  illustrate  the  significance  of 
the  freedom  to  choose  asymptotically  infinite  eigenvalues  and  eigenvectors 
to  achieve  desired  closed-loop  actuator  characteristics. 


Table  2.  F and  G Matrices 


-.  746 

. 387 

-12.9 

0. 

.952 

6.05 

.024 

-.  174 

4.31 

0. 

-1.  76 

-.  416 

. 006 

-. 9994 

-. 0578 

. 0369 

. 0092 

-. 0012 

F = 

1. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

-10. 

0. 

0. 

0. 

0. 

0. 

0. 

-5. 

0.  0. 

0.  0. 

0.  0. 

0.  0. 

20.  0. 

0.  10. 
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From  the  point  of  view  of  fighter  handling  qualities,  all  four  of  the  lateral 
cLxis  closed  loop  roots  have  desired  values  which  can  be  taken  from  MIL" 
F8785B,  as  done,  for  example,  in  reference  1 10],  The  results  are  as 
f ollows: 


(i) 


a)  Roll  subsidence  mode  -4.0 

b)  Dutch  roll  mode  -0.63  ± j 2.42 

c)  Spiral  mode  -0.05 

Each  of  these  poles  can  be  assigned  an  asymptotic  eigenvector,  x'"',  which 
distributes  the  modal  response,  e among  the  state  variables  and  outputs 
of  the  system.  However,  each  eigenvector  is  constrained  to  lie  in  a two 
dimensional  subspace  defined  by  Equation  (54).  An  element  in  this  sub- 
space was  selected  by  finding  the  best  linear  projection  of  an  unconstrained 
desired  vector,  x*,  on  the  subspace.  The  results  of  the  eigenvector 
selection  are; 

a)  Roll  subsidence  mode  (e 


Desired  x* 

= ll. 

0 

0 

1 

1 

V 

vj 

Attainable  x^^^ 

= ll. 

007 

0 

-.25  1 

.13 

-.  56J 

Dutch  roll  mode. 

real  part 

(e‘- 

63t^ 

cos  2. 

42t)x 

(2)^ 

Desired  x* 

= 10 

V 

1. 

0 j 

V 

vJ 

Attainable  x^^^ 

= lO  15.6 

1. 

0 1 

7.86 

-.103J 

c)  Dutch  roll  mode,  imaginary  part  (e  ‘ ^^^(sin  2.  42t)x^*^^) 


Desired  x* 


= [0  1. 


0 V vj 


Attainable  x 


= [0  1.  6.  16  0 -9.  49  14.  6j 


n c • 1 ^ / --05t  (4). 

d)  Spiral  mode  (e  x ) 

Desired  x*  = [v  v 0 1.  | v vJ 

Attainable  x^‘*^  = [-.  05  . 037  0 1.  0014  -.  0079  J 


A few  comments  are  in  order  to  explain  these  choices.  Consider,  for 
example,  the  roll  subsidence  mode.  The  desired  eigenvector  is  taken  to 
be  X*  = (1  0 0 V v)  which  means  that  the  mode  should  show  up 

dominantly  on  roll  rate,  but  not  on  yaw  rate  or  sideslip  (we  want  no  side- 
slip buildup  during  turn  entries).  These  are  basically  handling  quality 
considerations.  The  v's  in  the  vector  indicate  that  we  do  not  care  how 
much  of  the  mode  shows  up  on  these  components.  Certainly,  since 
0 = ^Pgdt,  some  mode  content  has  to  be  expected  on  x^  and  similarly,  if 
the  surfaces  are  actually  controlling  the  mode,  some  mode  content  should 

also  appear  in  x and  x . The  linear  projection  which  best  achieves 

^ ° (1) 

these  objectives  is  shown  as  x above.  Note  that  we  can  satisfy  our 
desires  almost  perfectly. 


Similar  arguments  also  apply  to  the  dutch  roll  mode.  In  this  case  we 
want  no  oscillatory  dutch  roll  content  on  roll  rate  and  bank  angle.  This  is 
a key  handling  quality  requirement  for  all  well-behaved  lateral  control 
laws  I 10 1 . 
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In  the  case  of  the  spiral  we  want  the  mode  to  show  up  dominantly  on  bank 
I angle  (corresponding  to  steady  turns)  with,  again,  no  substantial  sideslip 

component.  The  latter  is  a basic  turn  coordination  requirement. 

! 

Once  the  eigenvectors  are  specified,  it  remains  only  to  compute  the  matrix 
j via  Equation  (96)  and  then  to  select  eigenvectors  and  eigenvalue  ratios 

for  the  asymptotically  infinite  roots.  The  matrix  is  shown  in  Table  3. 

I 

Table  3.  C Matrix 
o 


-.  131 

. 567 

-.612 

. 160 

1.64 

-2.39 

. 0175 

. 0303 

1. 

0. 

0. 

1. 

For  decoupleci  actuators  the  matrix  N defined  in  (59)  and  matrix  S of  ec).  (60) 
were  chosen  to  be 


1 

o 

0 

1 

0 

N = 

0 

1 

s = 

0 

0.  5 

to  achieve  a two-to-one  bandwidth  ratio.  Then  from  (57)  and  (56), 


These  selections  complete  the  specification  of  the  performance  index 
except  for  the  scalar  pai'ameter,  c.  Optimal  controllers  were  computed 
for  several  values  of  this  parameter. 

Representative  results  are  shown  in  Table  4 and  Figure  2.  The  eigen- 
values and  eigenvector  characteristics  indeed  approach  their  predicted 
asymptotic  values,  as  can  be  seen  in  Table  4.  Reasonable  actuator 
bandwidths  are  obtained  with  p = 0.  0025. 

The  controller  gains  shown  in  Table  4 indicate  that  the  actuators  at  this 
design  value  are  essentially  uncoupled  (very  small  cross -feeds  from  one 
to  the  other)  and  have  just  enough  feedback  around  themselves  to 
produce  the  bandwidths  available  with  the  aircraft's  existing  hydraulics. 

Transient  response  characteristics  achieved  for  a range  of  p -values  are 
shov/n  in  Figure  2.  Although  only  roll  rate  and  sideslip  responses  are 
given  for  simplicity,  dramatic  improvement  in  cross  coupling  between 
roll  motion  and  dutch  roll  is  clearly  evident.  Roll  transients  approach 
their  intended  first  order  0.  25  second  response  time,  and  sideslip 
respon.ses  approach  tlieir  intended  sliglitly  oscillatory  second-order 
characteristics,  C.'ross  coupling  vanislies  as  the  eigenvectors  approach 
tlieir  intended  asymptotic  directions. 


This  example  illustrates  that  it  is  possible  to  satisfy  complex  handling 
quality  specifications  and  actuator  bandwidth  cliaracteristics  with 
standard  quatiratic  performance;  indic(;s  if  tin;  structure  of  the  criterion 
is  based  on  asymptotic  ctiaracf eristics. 


bt  ®*0  8 --  8f ’ * »ao 


Figure  2. 


'I'ransient  Responses 
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I'igure  2.  Transient  Responses  (concluded) 
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APPENDIX  A 


ZEROS  OF  MULTIVARIABLE  SYSTEMS 


This  Appendix  gives  a summary  of  the  important  concepts  and  an  outline 
of  the  literature  on  transmission  zeros  of  multivariable  systems.  While 
many  intuitive  notions  about  zeros  for  scalar  systems  are  retained,  the 
multivariable  case  is  not  a trivial  extension  of  the  scalar.  Examples  are 
given  to  illustrate  important  ideas. 

SUMMARY  OF  IMPORTANT  CONCEPTS 

1.  Given  the  system: 

X = Fx  + Gu 
y = Cx  + Du 

, „n  „m  p. 

(xeR.ueR  ,yeR) 

the  transmission  zeros  may  be  defined  to  be  the  set  of  complex 
numbers  \ which  satisfy  the  following  inequality: 


rank 


F -XI 
C 


G 

D 


< n f min  (m,  p) 


2.  Transmission  zeros  are  invariant  with  respect  to  coordinate  trans- 
formation, state  feedback,  and  nonsingular  input  and  output  trans- 
formation. 
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3.  Transmission  zeros  may  include  the  entire  complex  plane.  The 
remaining  points  will  assume  this  is  not  the  case. 

4.  Almost  all  nonsquare  systems  (m  k p)  have  no  zeros. 

5.  For  square  systems  (m  = p)  the  definition  in  l8j  is  equivalent  to 
defining  transmission  zeros  as  zeros  of: 

If  -XI  g1 


= det(F  - XI)  det[C  (XI-F)  G + Dj 

G,  For  square  systems,  the  number  of  zeros  is  at  most  n-m-d  where  d is 
the  rank  deficiency  of  CG. 

7.  Transmission  zeros  are  not  the  same  as  the  zeros  of  the  numerators 
of  the  various  transfer  functions  relating  inputs  to  outputs.  The  latter 
are  not  invariant  with  respect  to  feedback. 

8.  For  high-gain  feedback  systems  where  u = ky  the  system  poles  migrate 
to  the  transmission  zeros  or  infinity  as  k approaches  infinity. 


9.  A zero  at  X corresponds  to  the  system  blocking  the  transmission  of 
certain  input  signals  proportional  to  e^*'. 
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OUTLINE  OF  LITERATURE  ON  TRANSMISSION  ZEROS 


There  has  been  considerable  interest  recently  in  transmission  zeros  of 
multivariable  systems,  Rosenbrock  l8j  apparently  was  the  first  to  offer 
a definition.  Most  of  the  literature  offers  little  insight  and  many  of  the 
results  are  incorrect.  We  will  present  the  most  significant  results  here 
and  try  to  clarify  some  of  the  issues.  The  interested  reader  may  consult 
the  literature  for  details,  and  we  will  point  out  the  most  obvious  errors 
therein. 

Much  of  the  difficulty  arises  because  the  transmission  zeros  of  a system 
are,  in  general,  quite  different  from  the  zeros  of  the  various  scalar 
transfer  functions  relating  inputs  and  outputs,  as  in  the  single -input 
single-output  controllable-observable  case. 


Rosenbrock  [8j  defines  zeros  of  a transfer  function  matrix  in  terms  of 
the  numerator  of  the  diagonal  elements  of  its  McMillan  standard  form.  If 
y<s)  = G(s)  u(s),  (y  is  length  p,  u is  length  m,  and  G is  pxm),  then 
G(s)  = N(s)/d(s)  where  d(s)  is  the  monic  least  common  denominator  of  the 
elements  of  G.  Then  N(s)  can  be  brought  to  Smith  form  by  the  transforma- 
tion L(s)  N(s)  R(s)  = S(s)  where  L and  R are  unimodular  polynomial 
matrices.  Cancelling  any  common  factors  between  S(s)  and  d(s)  yields 
the  McMillan  form  of  G which  we  denote  by  M(s)  as  follows: 

M(S)  = ldiag(e  .(s)/»^/.  (s))  ! 0 J,  m >p  (A.l) 

I 1 I P»  ”P 


or 


or 


M(s)  = diag(e  ^(s)/i//.(s)).  m 


M(S)  = 


diag(e  ,(s)/l//^(s)) 


= p 

m < p 


0 

p - m , m 


(A.  2) 


(A.  3) 
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I 


The  poles  of  G are  the  zeros  of  the  denominator  polynomials  i/^^(s)  in  M(s) 
and  the  zeros  of  G are  the  zeros  of  the  e ^(s). 


The  following  illustrative  example  was  given  in  [8j. 


(s  + 1)' 


(s  + l)(s  + 2) 


G(s)  = 


(s  + l)(s  + 2) 


(s  + 2)‘ 


d(s)  = (s  + 1)^  (s  f 2)^ 


(s  + 2)‘ 


(s  + l)(s  + 2) 


N(s)  = 


(s  f l)(s  +2)  (s  f 1)  (s  4-  3) 


0 (s  + l)^(s  + 2)^ 


(s  + l)^(s  + 2)^ 


M(s)  = 


Notice  that  although  G is  proper,  M is  not.  Also,  G has  one  zero  at  -2, 
and  there  is  not  a zero  at  -3. 


I 
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Kwakernaak  and  Sivan  [7]  define  the  zeros  of  the  systena* 


i = Fx  + Gu  (A.  4) 

y = Cx  + Du 

(x  e u e y e R^) 

for  m = p and  D = 0 as  the  zeros  of  the  polynomial  det  (ls-F)x 
det(Cds-F)'^G).  IfH(s)  = C(sI-F)'^G  then: 

det[H(s)J  =-^.  d(s)  = det(sI-F)  (A.  5) 

Now: 

s"^  det[C(Is-F)'^G] 
det[Cs(sl-F)“^G] 

= det(CG)  <A.(i) 


lim  s’^^Ms)  _ lim 
js  1 - ® d(s)  ls|  - “ 

lim 


This  shows  that  the  degree  of  d(s)  is  greater  than  i/^s)  by  at  least  m,  hence 
i/Xs)  has  degree  n-m  (when  det(CG)  h 0)  or  less, 

Wolovich  [llj  attempts  to  present  a somewhat  more  general  state-space 
approach  to  multivariable  system  zeros  by  essentially  defining  a zero  as  a 
value  of  s for  which  the  matrix 


*ControUabUity  and  observability  will  be  assumed  throughout  unless 
otherwise  stated. 
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(A.  7) 


loses  rank.  The  normal  rank  of  a matrix  is  usually  taken  to  be  the 
maximum  rank  over  s and  a zero  can  be  thought  of  as  a value  of  s for  which 
the  rank  is  less  than  the  normal  rank.  Wolovich  incorrectly  states  that  if 
the  rank  of  G and  the  rank  of  C are  both  greater  than  or  equal  to 
r = min(m,p),  then  the  normal  rank  of  (A.  7)  is  n+r.  For  the  case  p = m 
and  D = 0,  Wolovich's  definition  is  equivalent  to  that  given  in  [7j  because 


= det(F-sI)  detlC(sI-F)  Gj 


(2.8) 


However,  it  is  possible,  as  the  following  example  demonstrates,  for 
observable -controllable  systems  with  C and  G each  of  full  rank  to  have  the 
matrix  in  (A.  7)  be  of  rank  less  than  n+r. 

Example;  Consider  the  controllable,  observable  linear  system  with  n = 3, 
m = p = 

X = fhc  + Gu 


where 


y = Cx 


0 1 0 

I-'  = 0 0 1 

0 0 ] 


G = 1 U 


1 0 u 

0 1 0 


They  point  out  that  a controllable  observable  system  with  G and  C full  rank 
may  have  zeros  everywhere  in  the  complex  plane.  They  call  these  systems 
degenerate  and  assume  that  all  their  systems  are  nondegenerate.  The 
above  definition  is  shown  to  coincide  with  those  of  Kwakernaak  and  Sivan 
and  Rosenbrock. 


It  is  also  noted  that 


rank 


= rank 


'f-xi  g' 

C D_ 

T(F  f GK  - xI)T 
U(C  + DK)T'^ 


TGV 

UDV 


(A.  11) 


for  any  nonsingular  T e nonsingular  V e nonsingular  U e R^^^ 

and  K e and  hence  that  the  transmission  zeros  {x]  are  invariant 

with  respect  to  coordinate  transformation,  state  feedback  and  nonsingular 
input  and  output  transformations. 


Davison  and  Wang  also  prove  that  the  class  of  square  systems  having  less 
than  n-m  zeros  must  lie  on  a hypersurface  in  the  parameter  space  for  those 
systems.  Thus,  they  conclude  that  for  almost  all  square  systems,  there 
are  n-m  zeros. 


Davison  and  Wang  have  many  other  control-related  results  on  zeros,  and 
claim  to  have  as  vvell  some  efficient  algoritlmis  for  calculating  zeros. 
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Desoer  and  Schulman  [15j  give  a very  interesting  account  of  zeros  using 

purely  algebraic  techniques.  The  rational  matrix  transfer  function  H(‘) 

is  viewed  as  a network  function  of  a multiport  and  H(s)  (n  xn.)  is  factored 
- 1 - 1 o 1 

into  D (s)  N.(s)  = N (s)  D (s)  . A zero  of  H(‘ ) is  defined  to  be  a point 

I I r r 

z where  the  local  rank  of  N^(z)  drops  below  the  normal  rank.  The  intuitive 

concept  that  a zero  of  a multiport  corresponds  to  the  system  blocking  the 

z t 

transmission  of  signals  proportional  to  e is  developed  rigorously  via 
several  theorems. 

MacFarlane  and  Karcanias  [16j  have  written  a rather  complete  survey  of 
poles  and  zeros  of  multivariable  systems.  They  treat  zeros  from  the 
algebraic,  geometric  and  complex -variable  points  of  view  and  give  several 
apparently  new  and  worthwhile  intuitive  approaches.  Furthermore,  there 
are  no  glaring  errors  in  the  paper.  This  is  definitely  good  (though  lengthy) 
reading  on  the  subject. 
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